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We analyze the signal-to- noise ratio for a relic background of scalar gravitational radiation composed 
of massive, non-relativistic particles, interacting with the monopole mode of two resonant spherical 
detectors. We find that the possible signal is enhanced with respect to the differential mode of the 
interferometric detectors. This enhancement is due to: (a) the absence of the signal suppression, for 
non-relativistic scalars, with respect to a background of massless particles, and (b) for flat enough 
spectra, a growth of the signal with the observation time faster than for a massless stochastic 
background. 
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The detection of a background of relic gravitational ra- 
diation is undoubtedly one of the main challenges of the 
physics of the XXI century. The detection of, or even a 
direct experimental bound on, a stochastic gravitational 
background of primordial origin would considerably im- 
prove our knowledge of the very early state of the Uni- 
verse and our understanding of physical processes at a 
Planckian energy scale (see for instance jjj an references 
therein) . 

In view of the present and forthcoming data, soon 
available from the cross-correlation of all existing (res- 
onant and interferometric) gravitational antennas, the 
signal-to-noise ratio (SNR) due to a stochastic back- 
ground has been analyzed in full details (see 0] an ref- 
erences therein), in particular for a background of (ul- 
tra)relativistic particles, whose spectral distribution in 
frequency is identical to their momentum distribution. 

As far as the tensor sector of gravitational radiation is 
concerned, the relativistic assumption is certainly appro- 
priate and self-consistent. However, it could not be fully 
appropriate in order to analyze the contribution of a pos- 
sible scalar sector, which could correspond to a stochas- 
tic background of massive particles (for instance, dilatons 
pj), possibly non-relativistic at the time of present ob- 
servations. 

A non-relativistic background of massive scalar waves 
requires indeed a generalization of the standard SNR ex- 
pression which takes into account the relevant mass pa- 
rameter of the energy spectrum in momentum space [0 . 
In that case it has been recently shown that, in spite of 
the induced suppression factor (with respect to massless 



particles) (p/E) 4 <^C 1, it is possible to obtain a resonant 
response also from the non-relativistic part of the scalar 
wave spectrum, provided the mass lies within the sen- 
sitivity band of the detectors []I|g. Also, if the energy 
density of the non-relativistic scalar background is a sig- 
nificant fraction of the critical energy density, the signal 
could be detectable by future planned detector configu- 
rations, such as "Advanced LIGO" 0. 

The analysis of plE) refers however to the SNR pro- 
duced by a massive scalar background in the cross- 
correlation of two interferometric detectors. In partic- 
ular, for a non-relativistic scalar wave, the suppression 
factor which appears in the pattern function is due to the 
interaction of the scalar wave with the differential mode 
of the interferometers, which is characterized by a 3 x 3 
traceless response tensor. As already pointed out in ||, 
such a suppression may be absent in the monopole mode 
of resonant spherical detectors, whose response tensor is 
proportional to the 3x3 identity matrix M . In addition, 
in the resonant response of a sphere to a non-relativistic 
scalar background, and for flat enough spectra, the in- 
frared cut-off of the SNR integral is provided not by the 
noise power spectrum of the detector (as usual in the case 
massless particles), but by a minimal frequency scale de- 
termined by the mass of the background field and by the 
total observation time, T. As a consequence, the corre- 
sponding SNR may rise with T much faster than usual 
(depending on the power of the scalar spectrum) . 

The aim of this paper is to present, and briefly discuss, 
these two effects (i.e., the absence of non-relativistic sup- 
pression, and the unconventional T-dependence of SNR) 



which are new (to the best of our knowledge), and which 
could greatly favour the detection of a non-relativistic 
scalar background through the correlation of two (or 
more) resonant spheres. 

We start by recalling the general expression for the 
optimized SNR obtained through the cross-correlation 
of two detetctors, with one-sided noise spectral den- 
sity Pi(|/|), P2G/I), an d generated by a stochastic back- 
ground of scalar particles with mass m and spectral en- 
ergy density fl(p) in momentum space [Mpl: 
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Here T is the total observation time, rh = m/27r, and Hq 
is the present value of the Hubble parameter, which ap- 
pears in the above equation because the energy spectrum 
p{p) of the background is expressed in units of critical en- 
ergy density Hif: 
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where Mp is the Planck mass. Finally, in eq. (Q), j(jp) 
is the overlap reduction function [yjj^], which accounts 
for the relative orientation and location of the two grav- 
itational antennas, and determines their efficiency in the 
detection of a given background. For a plane wave of 
momentum p = rip, where n is a unit vector specifying 
the propagation direction, 7(7?) is defined by 
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where the integration is extended over the full solid angle 
and the normalisation constant has been chosen so that 
- in the massless case - one obtains j(p) = 1 for the 
differential mode of two coincident and coaligned inter- 
ferometers. Here x\,X2 are the positions of the centers of 
mass of the two detectors, and F l ,F 2 are the so-called 
antenna pattern functions pLQ|, determined by the wave 
polarization tensor, e ab (n) (which depends on the spin- 
content of the background) , and by the detector response 
tensor D ab (which depends on the geometrical shape of 
the detector). In general: 



F l = q t e ab (n)D ab 



(4) 



where % is the effective coupling of the background to 
the detector, normalized to 1 for the geodesic coupling 
to metric fluctuations MM- 

Note that, for a better comparison with experimental 
observables, we are choosing "unconventional" units h = 
1, so that energy and frequency simply coincide, E = 



(p 2 + rh 2 ) 1 / 2 — /. Note also that, for m = 0, p = / and 
one recovers the standard form of SNR appropriate to a 
background of massless particles pl|| ; the usual result for 
a background of relic gravitons then follows by inserting 
into eqs. (ra), (H) the appropriate spin- two polarization 
tensor, with q = 1. 

For a massive scalar wave, however, one should in gen- 
eral distinguish two different pattern functions, one re- 
lated to the geodesic coupling of the detector to the scalar 
component of metric oscillations M, the other related 
to the direct, non-geodesic coupling of the scalar charge 
qi of the detector to the spatial gradients of the scalar 
background |9| . By using the transverse and longitudinal 
decomposition of the polarization tensor of the massive 
scalar wave with respect to the propagation direction n, 
and defining 



T a b = (Sab- n a fib), L ab = h a h b , 



(5) 



the two patter functions can be written, respectively, as 
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The first one applies when considering the response of 
the detetctor to the spectrum of scalar metric fluctua- 
tions induced by the scalar background, the second one 
when considering the direct response of the detector to 
the fluctuations of the background itself. For the differ- 
ential mode of an interferometric antenna, in particular, 
D ab S a b = and, in both cases, the pattern function of 
a massive scalar wave turns out to be proportional to 
the massless one, but with the strong suppression factor 
(p/E) 2 for non-relativistic modes ||. 

A resonant sphere, on the contrary, has a monopole 
mode characterized by the response tensor D ab — S ab [J7) , 
so that D ab L ab = 1, D ab T ab = 2. The non- geodesic part 
of the corresponding pattern function, F q (n), is still sup- 
pressed for all non-relativistic modes. For the geodesic 
part, however, the angular dependence completely disap- 
pears, 
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and the overlap function (|3|), for the monopole modes of 
two spheres, takes the simple form 
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where d is the spatial distance between the detectors. 
The response to non-relativistic modes {p <C m) is no 
longer suppressed (actually, their response is enhanced by 
the factor 9/4). This is a first interesting result, which 



already at this level shows that a resonant sphere is a 
promising device for the detection of a non-relativistic 
background of massive scalar particles. 

The second interesting result follows from the study of 
the integral (Q) for the SNR associated to the monopole 
modes of two spheres, and induced by a cosmic back- 
ground whose spectrum f2(p) is dominated by the non- 
relativistic sector p < m. Such a relic scalar spectrum 
is naturally expected in string cosmology models of the 
early Universe |L0[ and, as discussed also in M , it can be 
simply approximated by the power-law behaviour 
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where f^o is the present fraction of critical energy density 
stored in the scalar background. The low energy branch 
of the spectrum is assumed to be non-decreasing (S > 0), 
to guarantee a finite energy density, 



dhip £l(p) = f^Q. 



(10) 



More complicated distributions are also possible, with 
O(p) peaked at a scale p m < rh (see plfLOfl), depending 
on the details of the cosmological mechanism of produc- 
tion. But the simple monotonic spectrum (pj) is already 
appropriate to the illustrative purpose of this paper, and 
will be adopted in what follows for all our estimates of 
the sensitivity. 

Inserting into eq. (|lj) the expression for the overlap 
reduction function given by eq. (|j) and the spectrum (g) , 
and assuming (for a resonant response) that to lies within 
the detector sensitivity band, so that Pi(rh) are finite, the 
resulting SNR would seem infrared divergent for 5 < 1. 
Indeed, for p — * 0, ~/(p) — * const, Pi — » const, and 
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Unlike for massless particles, the infrared part of the inte- 
gral is not killed by the instrumental noises, whose power 
spectra Pi(\f\) (appearing in the denominator of eq. (Q)) 
blow up to infinity when their argument approaches zero. 
In our case, when p — ► 0, Pi keep frozen at the frequency 
scale fixed by the mass of the background. 

It must be noticed, however, that the presence of a 
lower bound fixed by to in the allowed range of frequen- 
cies, / = (p 2 + TO 2 ) 1 / 2 > to, toghether with the existence 
of a "maximal" finite time scale, determined by the obser- 
vation time T, necesssarily implies an intrinsic infrared 
cut-off for the momentum variable in the SNR integral. 
Indeed, for a given T, the minum observable frequency 
interval 



Af=(p 2 +m 2 ) 1/2 -m>T-\ 
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defines, in the limit p — > 0, the minimum momentum 
scale 



Pn 



= {2m/Tf 2 , 
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to be used as the effective lower bound in the SNR in- 
tegral. This instrinsic cut-off may thus introduce an im- 
portant T-dependence in the final value of SNR. 

In order to illustrate this effect, and to provide some es- 
timates, let us assume for simplicity that the instrumen- 
tal noises P\ and Pi are equal, and nearly constant, for p 
ranging from to to (i.e., for frequencies between to and 
V2to): Pi (to) = P-z{rn) = Sh(rh). Also, we shall assume 
that the distance between the two spheres is negligible. 
The SNR, for a background of scalar metric fluctuations 
charaterized by the spectral distribution (|9|), interacting 
with the monopole modes of two spheres, then reads 
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For 5 < 1 the integral is dominated by the infrared cut- 
off Pmin (see eq. (fill)). We obtain the following estimates: 
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For 5 > 1 the integral is instead dominated by the end 
point p — rh of the spectrum, and the SNR reads 
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is a number of order one. The three different cases give 
results which are very similar, except for the dependence 
on the observation time T, and it is now evident that, 
for a given T, the sensitivity is strongly enhanced for flat 
enough spectra (5 < 1). 

It is important to stress that such a dependence of the 
SNR on the total observation time yields an effective en- 
hancement - with respect to the usual square root law - of 
the signal, for a massive stochastic background, only if to 
is not much smaller than the typical resonant frequency 
/o of the detector. Indeed, if to <C /o, the infrared cut- 
off (O) is ineffective, being suppressed by a very high 



instrumental noise. It is also worth stressing that this 
enhancement is present only in the correlation of two 
spheres. Indeed, only because of their particular overlap 
function (||) the SNR integral may be dominated by the 
infrared cut-off, and not by the peak of the spectrum (as 
in the case of interferometers) . 

For a numerical estimate of the sensitivity we shall 
now compute the fraction of critical energy density Qo 
required to have a detectable signal, say SNR > 5. 
We take T ~ 10 7 sec as a typical obervation time, and 
Sh ~ 10~ 46 Hz as a possible realistic noise density for 
a resonant sphere fll]Jl3] at a frequency / = 3 x 10 3 Hz. 
The minimal required energy density, for 8 = 1/2, 6 = 1 
and 8 > 1 is given, respectively, by 
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suppression for the monopole mode of the spheres and, 
for a flat enough spectrum, the faster increase of the sen- 
sitivity with the observation time. The numerical esti- 
mates reported in this paper refer to the case in which 
the mass of the scalar particle is within the bandwidth 
of the antennas. The bandwidth depends on the noise 
performance of the readout system, and is limited to a 
few Hz in the currently operating resonant bars pjfl . By 
employing a xylophone of advanced hollow spheres p2| , 
or the recently proposed dual spheres p3, it should be 
possible however to open up the bandwidth and scan a 
wide mass range at a high level of sensitivity, compara- 
ble with the present numerical estimate, thus providing 
new important information on possible ultralight scalar 
partners of the gravitons, and on fundamental string/M- 
thcory models of the early Universe. 
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where we have used Ho = 3.2 ho x 10~ 18 Hz. 

We may note, for comparison, that eq. (|2l|), with to 
replaced by the typical resonant frequency /o, and F(8) 
always of order one (but replaced by a different func- 
tion of 8) also describes the sensitivity of two spheres 
to a massless stochastic background, for all values of the 
spectral index. Eq. (|2l|), with a different, slightly smaller 
F(8), also estimates the sensitivity of two coincident and 
coaligned interferometers to the massive spectrum (B), 
for any 8 (even if the maximal sensitivity is expected 
at frequencies lower than 1 kHz). For the interferome- 
ters, however, the sensitivity is to be further depressed 
if the spectrum is peaked at scales smaller than to, be- 
cause of the non-relativistic suppression factor B. For 
the spherical detectors such a suppression is absent, and 
the sensitivity remains at the above levels independently 
from the shape of the spectrum. 

We may therefore conclude that the cross-correlation 
of two resonant spherical detectors seems to be, in princi- 
ple, particularly appropriate to the search for a possible 
relic cosmic background of non-reativistic scalar parti- 
cles, because of two effects: the absence of non-relativistic 
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